Fluctuations and Tc reduction in cuprate superconductors 
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We analyse fluctuations about Tc in the specific heat of (Y,Ca)Ba2Cu307_i, YBa2Cu408 and 
Bi2Sr2CaCu2 0g+5 . The mean-field transition temperature, TJ"^, in the absence of fluctuations lies 
well above Tc especially at low doping where it reaches as high as 150K. We show that phase and 
amplitude fluctuations set in simultaneously and T™^ scales with the gap, Aq, such that 2Ao/kBT^-^ 
is comparable to the BCS weak-coupling value, 4.3, for d-wave superconductivity. We also show 
that is unrelated to the pseudogap temperature, T*. 

PACS numbers: 74.25.Bt, 74.72.-h, 74.62.Dh, , 74.62. Fj 



Many authors have suggested that pairing in the 
cuprates begins well above Tc- Emery and Kivelsen 



argued that the low superfluid density. 



the 



cuprates leads to phase fluctuations below the mean-field 
(MF) transition temperature, T™^, resulting in a phase- 
incoherent state with a finite pairing amplitude [l|. In 
this view phase coherence is not established until a lower 
temperature, the observed T^. Experimental support for 
this picture may be found in the high-frequency optical 
studies by Corson et al. 0] . Additionally, the underdoped 
normal state (NS) exhibits a depletion of the density of 
states (DOS) near the Brillouin zone boundary due to 
the presence of a pseudog api. The pseudogap seems to 
close abruptly at Pcrit=0.19 holes/CuQ. Some authors 
have drawn these two strands together proposing that 
the pseudogap corresponds to thephase-incoherent pair- 
ing state between and TJ"/[i[|]. The pseudogap T* 
line (below which pseudogap effects are observed) would 
then correspond to the doping-dependent T™/. 

Here we analyze the fluctuations in specific heat, Cp, 
to calculate T™-^ and the mean-field jump, A'j"''^ , in spe- 
cific heat coefficient, 7. We find that at all doping levels 
T^""^ lies well above the observed Tc, reaching as high as 
113K for (Y,Ca)Ba2Cu307_5 (Y,Ca-123) and 150K for 
Bi2Sr2CaCu2084_5 (Bi-2212). Our approach is similar to 
that of Meingast et al. \^ using thermal expansion data. 
But, where they identified the pair- fluctuating state with 
the pseudogap, we show they are distinct. Further, we 
flnd that T™-^ and T* are very different in magnitude and 
doping dependence, conflrming that they are unrelated. 
Over a wide doping range T™^ scales with the supercon- 
ducting gap such that 2Ao/kBT™f remains comparable 
to the weak-coupling BCS d-wave value of 4.3. 

Even without such an analysis the idea that the pseu- 
dogap is a phase incoherent pairing state faces an insur- 
mountable obstacle. If the pseudogap arises merely from 
thermal phase fluctuations then at T=0 there should be 
no remnant pseudogap effects. But, even at T=0 the 
pseudogap weakens the SC ground state, abruptly reduc- 
ing the condensation energyfil and superfluid density [tI, 



as doping is reduced below Pcrit- The pseudogap thus 
coexists with SC at r==0 and must be distinct from fluc- 
tuation effects above Tc. 

Firstly, we show that phase and amplitude fluctuations 
set in simultaneously. Emery and Kivelsen[l| deduce that 
phase fluctuations become important above a tempera- 
ture Tg given by fc^Tg ~ AVq where A ^ 1 and Vq is the 
phase stiffness, Vq = a?i^ns(0)/4m*. The length scale, a, 
was defined as a = -^Tr^ for isotropic 3D behavior and 
a — max((i, \/7r^j^) for 2D where d is the mean interlayer 
spacing. Vq is related to the penetration depth, Xab, viz: 
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The condensation energy, Uq, is given by 
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where the last equality comes from Ginzburg-Landau the- 
ory. Combining these, we find 



keTe - AVo = [AA/ti) Uo^{Q), 



(3) 



where £7(0) — Tr^^^a is the coherence volume for Cooper 
pairs. Following BulaevskiiQ we adopt the criterion for 
amplitude fluctuations as 

keTa^p ~ Uo^iQ), (4) 
which leads immediately from Eq. [3] to the relation 



keTe ~ AVa - (4A/^) Ua^{Q) 



{AA/t:) ksTamp- 



(5) 



As A « 0.9 for 2D[i| then the conditions for phase and 
amplitude fluctuations are equally restrictive. For a ho- 
mogeneous system they both must set in simultaneously. 
We thus question the widely-accepted phase-fluctuation 
model of Emery and Kivelsen[l| and its implementation 
by Corson et al.^. If Tg and Tamp greatly exceed T™^ 
then the transition occurs essentially at T^f . But, if Tg 
and Tamp are comparable to or less than T™^ (as is the 



2 



case) then will be suppressed below T™^. Between Tc 
and TJJ^f both amplitude and phase will fluctuate, not 
just the phase. It is our aim to determine how large this 
Tc suppression is. 

The fluctuations in C'p{T) have been analyzed[^ 
by separating Cp into a symmetric fluctuation term, Cp , 
and an asymmetric MF term, C^^'^. In the 3D-XY 
model Cp near Tc may be approximated by 




(t = (T/Tc 
(t > 0). 



is the MF step at Tc and A' 



1)<0) 



(6) 

A+ = 



4A:B/(97r2f7(0)) While Eq. [6] is not strictly correct 

deep in the critical region it captures all the main phys- 
ical features of the more complex crossover from critical 
to MF behavior [l^ and e.g. accurately represents the 
critical behavior of He^ at the superfluid transition ^9]. A 
plot of ACp versus In |i | gives two parallel lines offset 
by AC™-'^'°. In practice this plot exhibits negative cur- 
vature for sufflciently small \t \ due to minor transition 
broadening. The effect of the resulting spread in Tc was 
modeled by replacing thy t* = [t^ + At^)^/^ in the above 
expressions for ACp [l3l |. 

For Bi-2212 ACp was found to collapse rapidly with 
the opening of the pseudogap at Pcrit-, falling to zero 
near optimal doping p = 0.16 holes/Cu[lOj. Below this, 
Cp{T) is dominated by fluctuations alone and is symmet- 
rical about Tc- This is puzzling because the specific heat 
jump should remain finite, consistent with the second- 
order phase transition. We resolve this below. 

Fig. [ija) shows ^{T) reported by Loram et al. Q] for 
Yo.8Cao.2Ba2Cu307-5 with (5=0.25 and p=0.186.' The 
dashed line is ^n{T)- Because entropy S — J 'ydT there 
are two entropy balance conditions: (i) the area abc 
equals the area cde. This helps to establish the T- 
dependence of 7„ below Tc- In this case there is no pseu- 
dogap and jniT) is constant across the entire T-range. 
For lower doping where the pseudogap is present we use a 
triangular gap which fills with increasing temperature [Sj]: 

7n(T) = 7„(oo) [1 - tanh(?9) In (cosh(i9))] , (7) 

where = Eg/ikpT . The second entropy balance con- 
dition concerns the fluctuation term which reduces Tc 
below TJ^f . Thus, (ii) the entropy equal to the forward 
cross-hatched area between Tc and T™^ equals the fluctu- 
ation entropy given by the backward cross-hatched area 
under the fluctuation term, 7^' , which includes both crit- 
ical and Gaussian fluctuations. That is 



{iT^ - In) dT, (8) 



This construction enables TJ"-^ to be estimated. Further- 
more, the apparent mean- field step, Aj™f''^, at Tc is also 




_i I I L 



20 40 60 80 100 120 



temperature [K] 

FIG. 1: (Color online) (a) Analysis of -/(T) for 
Yo.8Cao.2Ba2Cu3O6.75 showing the MF component, ■yT^ , and 
the fluctuation component above and below Tc- Deduced val- 
ues are rc=82.79K, rr-^=97.50K, A7'°'=3.00, A7™-'^=1.51 
and A7'"^'''=2.31 mJ/g.at.K^. (b) a similar analysis for 
YBa2Cu408 showing 7™-'^ (blue curve) and the pseudogapped 
7„(r) (red dashed curve). The two additional 7(r) curves for 
2% Zn (cyan) and 4% Zn (green) coincide with 7„ (T) and thus 
confirm our pseudogap model. The upturns at low T are due 
to impurities. Deduced values are Tc^Sl.OOK, rj"^=91.92K, 
A7'°*=l.ll, A7'"^=0.90 and A7'"-'^^''=0.52 mJ/g.at.K^ 



smaller than the (hypothetical) mean-field step, A7™-'', 
that would occur at TJ"^. These jumps are defined in 
the figure. (The same superscript notation is used for 
the jumps in specific heat ACp*, ACp''^'*^ and ACp''^). 

We proceed as follows. We combine the first entropy 
condition with Eq. [7] to establish 7„ (T) . We then plot 



Cp above and below Tc using Eq.[6]to determine AC™ ■^'^ 
and hence A7™-^'° [l3| . We then construct a power-law fit 
to 7s (2^) at low T that reproduces this value of A7™^'°. 
This is j™-^{T) which is slightly superlinear consistent 
with the predominant d-wave gap structure. Lastly, we 
impose the second entropy condition (Eq. [5]) to deduce 
T™-^ and A7™-''. Naturally there are errors inherent in 
such a construction. But while they grow with under- 
doping they do not impact on any of our conclusions. 

In the example shown in Fig.[l]rc is reduced by (14. 7± 
1.5) K below T^f. At the same time A7™/ = (2.31 ± 
0.08) mJ/g.at.K^, significantly more than A7'"-^'°=1.51 
mJ/g.at.K^. The analysis was carried out for many sam- 
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FIG. 2: (Color online) The doping dependence of evaluated parameters, (a) and (b) show Tc, T"'^, A7™-^ and A7'"^'° for 
Yo.8Cao.2Ba2Cu307_i ; (c) and (d) show the same for Bi2Sr2CaCu208+i . In (a) values are also shown for YBa2Cu307_i 
(arrows). The SC gap, Ao (red open triangles), and pseudogap Eg (blue crosses) from ref. 16) are also plotted, scaled by the 
factor (l/2.5fc_g). (c) also shows Ao values for Bi-2212 from B\g Raman (red open triangles), (e) shows Te, T"^^ and T*[2l[. 



pies with different oxygen contents. Values of T^. and 
T™-^ are plotted in Fig. [2ja) along with values of 
and A7'"^'° in panel (b). A similar analysis was done on 
specific heat data for Bi-2212[6'|. Here, instead of using 
Eq. [7] the full bilayer ARPES dispersion was usedTJJ, 
thus incorporating the van Hove singularity and pseudo- 
gap. The pseudogap was implemented as beforejl4| using 
a finite- Fermi-arc model. Results are plotted in Fig.l^^c) 
and (d). Several key conclusions can be made: 

(i) Like Meingast et al. jsj] we find T™^ (p) continues to 
rise with decreasing doping and only falls at the lowest 
doping levels. Underdoped samples show a reduction in 
Tc below T^f as large as 35-40K for (Y,Ca)-123 and 60K 
for Bi-2212, reflecting the larger anisotropy in the latter 
compound. The shift is also large for pure YBa2Cu306.97 
with Tc = 92.9K and T^^ = 112.3K (see arrowed data 
points Fig. [DJa)). 

(ii) while A7™-''''' = at lower doping (and the specific 
heat anomaly then becomes a pure symmetric fluctuation 
term) A7™-^ remains finite in the absence of fluctuations 
and may only reach zero near the onset of SC at p « 0.05. 
This removes the puzzle of the seemingly zero MF step. 

(iii) Panels (a) and (c) show the pseudogap energy. Eg , 
as previously reported[g,[l4|. Coincident with the abrupt 
opening of the pseudogap at p « 0.19 there is an abrupt 
reduction in all values of A7 showing that, even after 
removing fluctuation effects, the pseudogap still plays a 
decisive role in weakening the condensate and reducing 
the density of Cooper pairs. 

(iv) We compare T™^ with the r=0 SC gap, Aq, for 
the two systems in Fig. [2][a) and (c), respectively. For 
(Y,Ca)-123 values of Aq are from the speciflc heat[6| and 
for Bi-2212 from the Raman Big gap[l5|. In both cases 



2Ao/kBT^^ « 5 across the entire overdoped region, lit- 
tle more than the d-wave mean-field BCS value of 4.3. 
The old puzzle that 2Ao/fcsTc increases steadily with de- 
creasing doping[15] is now resolved by referencing to T™-^ 
rather than Tc- Eventually, at low doping T™f{p) falls 
below Ao/2.5kB due to the pseudogap progressively re- 
moving spectral weight. If it were not for the pseudogap 
T'P^ {p) would probably track Ao/2.5fcB across the entire 
SC domain. The fact that A7 is immediately reduced at 
p = 0.19 while TJJ^^ is only gradually reduced by the 
pseudogap is precisely what is expected in a competing 
gap scenario [l6|. Thus Tc is reduced both by fluctuations 
and eventually by the pseudogap once Eg becomes com- 
parable to Aq . A corollary is that Tc,max , the maximum 
in the Tc{p) phase curve, has no physical significance. It 
is the r=0 rf-wave SC gap, Ao, which is the truly fun- 
damental quantity and in a BCS scenario TJ"-^ will scale 
with Aq until the pseudogap opens, as observed. 

(v) If we use An values for (Y,Ca)-123 from infra-red c- 
axis conductivity]!!] we obtain 2Ao/A:bTJ"/ « 4.2 - 4.4 
in even better agreement with the BCS weak-coupling 
value. In fact, it is probably fair to state that gap mag- 
nitudes are not known sufficiently accurately to discount 
precise agreement with the weak-coupling value. 

Fig. [ijb) shows a similar fluctuation analysis on un- 
published 7(r) data for YBa2(Cu,Zn)408 with 0, 2 and 
4% Zn on the planar Cu sites. This shows the rapid 
suppression of both Tc and the jump in speciflc heat 
due to impurity scattering. The high rate of suppression 
dTc/dx — 13 K/%Zn is typical of underdoped cuprates 
and reflects the presence of the pseudogap which en- 
hances the pairbreaking scattering rate due to the re- 
duced DOS 18']. We use Eq. [7]to flt the pseudogapped 
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7n {T) (shown by the dashed red curve) and the fit is con- 
firmed by the 2% and 4% curves for 7(T) in Fig.[T]^b) for 
which the NS values coincide with the dashed red curve. 
The upturns in "f(T) at low T are due to a small fraction 
of impurity and need not concern us. 

Next, the values of 7^-^(T) (blue curve) are determined 
by fitting a power law to 7s(T). The complication of 
the upturn in the experimental data at low T is averted 
by insisting on entropy balance such that the area be- 
tween the dashed red curve and blue curve below the 
crossing temperature, Tcross = 61. 2K, equals the area 
between the black and dashed red curves above Tcross- 
We thus obtain T^f = 91.92K from = 81.00K; while 
A7*°*=l.ll, A7"-'^=0.90 and A7"-'^'0=0.52 mJ/g.at.K^. 
The depression in Tc due to fluctuations is ATc = 10.92K, 
rather less than the value ATc — 33. 5K obtained for 
Yo.8Cao.2Ba2Cu307_5 at the same doping state. This 
shortfall is probably due to the enhanced superfluid den- 
sity in YBa2Cu408[7| which, according to Eqs. [T] and 
[21 will suppress fluctuations. This implies that the gap 
magnitude is less in the latter compound, perhaps due 
to the proximity effect between CU2O2 chains and Cu02 
planes which is expected to lower the SC gap magnitude. 

There is little data available on the gap magnitude in 
YBa2Cu408 but Janossy et a/. [13] have carried out pre- 
cise measurements of the T-dependence of the spin sus- 
ceptibility below Tc using Gd electron spin resonance. 
They find an excellent MF d-wave fit with maximum 
(antinodal) gap Ao = 190K, giving 2Ao/fc_Brc=4.75. How- 
ever, by referencing to T^^ we obtain 2Ao/fcsT™-^=4.14, 
now very close to the weak-coupling value. 

Finally, Fig. [DJe) compares the various relevant tem- 
perature scales, Tc{p), T"'-f{p) and the pseudogap line 
T*{p) for Yo.8Cao.2Ba2Cu307_5. A central conclusion 
is that, contrary to some authors [l^] the pseudogap line 
T*{p) does not merge on the overdoped side with Tc{p), 
still less with the more fundamental quantity T^f{p). 
This question was effectively put to rest several years 
ago 21j by an extensive study of the resistivity of high- 
quality thin films of (Y,Ca)Ba2(Cu,Zn)307_5. The com- 
bination of Zn substitution and high magnetic fields al- 
lowed Tc to be suppressed so as to expose the evolu- 
tion of T* below Tc- Panel (e) reproduces these values 
of T*{p) (blue data points; solid=films; open=sintered). 
They extend below the unperturbed Tc value, descend- 
ing towards zero at p « 0.19. The solid blue curve is 
a power-law fit consistent with a terminating quantum 
critical point 22] ■ Importantly Zn substitution and mod- 
erate magnetic fields do not modify T* [2T] while they do 
suppress Tc- In this way it is straightforward to distin- 
guish between pseudogap effects and SC fluctuation ef- 
fects in the transport properties. T* evidently dips well 
below T™^ . The temperature scales shown in Fig. [2fe) 
are all of comparable magnitude and it is therefore not 
surprising that they have been confused in the past. 



In summary, we have carried out a fluctuation analysis 
of specific heat data to determine the mean-field transi- 
tion temperature, T™^, and the mean-field jump in spe- 
cific heat coefficient, A7™-^. T^^ rises rapidly above Tc 
with decreasing doping, reaching values of about llOK for 
YBa2Cu307_5 and Yo.8Cao.2Ba2Cu307_5, and as high 
as 150K for Bi2Sr2CaCu208+5. This shows the funda- 
mental importance of fluctuations in the cuprates. A7™^ 
remains non-zero across the phase diagram, as it must for 
a second-order phase transition. The long-standing puz- 
zle that 2Ao/fcBTc grows with reducing doping is resolved 
by replacing Tc by TJJ^^ . Across much of the SC phase 
diagram 2Ao/fcsT™-'' remains close to the weak coupling 
BCS value. T* is shown to be distinct from T^-f . 
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